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ABSTRACT 


The first-stage, powered-flight trajectory of a large rocket pow- 
ered vehicle is studied by varying the initial acceleration, the vertical 
fiers time, and the initial tilt angle. Trajectories were computed on an 
IBM 650 digital computer. Specific areas of interest with respect to high 
initial acceleration rockets are the feasibility of using the "gravity turn" 
maneuver to obtain low burnout flight path angles, and the determination of 
maximum energy trajectories for various values of intial acceleration. 


mesults indicate that a relatively low initial tilt angle followed 
by a "gravity turn" maneuver is not adequate to achieve low burnout flight 
path angles for high initial acceleration vehicles. For values of initial 
acceleration of about 2.5 to 3.0 a large percentage of burning time is spent 
in the programmed tilting phase, which results in lift load factors of the 
250? of .Ö to 1.2. 


Maximum energy trajectories occur at specific values of burnout 
flight path angle for the initial accelerations considered. These burnout 
angles start at about fifty-five degrees for an initial acceleration of 3.0 
and decrease to approximately zero degrees for an initial acceleration of 1.5. 


Burnout conditions of velocity, altitude, energy, and flight path 
angle are plotted for the trajectories computed. The trajectories most closely 
approximating the maximum energy cases are included in tabular form. 


Thesis Supervisor: Paul E. Sandorff 


Title: Professor of Aeronautics and 
Astronautics 
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OBJECT 


The object of this thesis is to study the early powered-flight 
trajectory of a large rocket powered vehicle. The effects on the first- 
stage trajectory of varying vertical flight time, initial tilt angle, 
and initial acceleration, are of primary interest, especially as they 
affect the maximum burnout energy conditions. 

Of interest also is the feasibility of using a relatively small 
initial tilt angle followed by a "gravity turn" to reach practical burnout 
conditions of velocity, altitude, and flight path angles for vehicles with 


nign initial acceleration. 





CHAPTER 1 


INTRODUCTION 


The study reported herein is concerned primarily with the initial 
portions of the powered-flight trajectory of a large, single stage, rocket 
powered vehicle. Conceptually, this vehicle could be the booster stage of 
an ICBM or a satellite launcher. 

Usually there are three phases to the initial flight trajectory 
of a large ballistic missile or satellite launching vehicle. These phases 
include a vertical flight phase, a tilt phase, and a gravity turn phase. 

The gravity turn phase is customarily followed by a period of "constant- 
attitude thrust", during which the major portion of the flight velocity is 
5277 66 (his latter regime is not considered in this study, 

A vertical launch for a large, rocket-povered vehicle of current 
design is necessary due to the inability of the vehicle structurally to vith- 
stand the transverse loads which would be present during an inclined launch. 
Nerbical or near vertical flight is also necessary in order to achieve 
2015: 006. Usually the vertical flight path is followed for a short time, but 
che time of vertical flight must be carefully selected in order to achieve a 
trajectory which minimizes propellant expenditure. 

Upon completion of the vertical flight phase, a tilting phase is 
commenced. Tilting is normally accomplished by deflecting the thrust vector 


Of the vehicle to produce a tilting moment according to some selected program; 





this changes the attitude of the vehicle, and subsequent thrusting changes 
the velocity vector. This maneuver is non-optimum and is best completed 
quickly; however, the tilt rate must not be so rapid as to exceed practical 
Limitations of the vehicle control and structure. The tilting phase is com- 
pleted when the vehicle body axis and the thrust vector are both aligned 
with the vehicle velocity vector. 

The third phase of the conventional trajectory concerns the flight 
regime where this alignment exists and a relatively slow turning path follows, 
brought about by the component of gravity transverse to the flight path. 

This phase hopefully terminates at an altitude above the sensible atmosphere, 
and with an attitude that matches the subsequent constant-attitude thrust re- 
gime in such a manner that the best overall trajectory performance is obtained. 

The important problem of proceeding from the earth's surface, 
through the three phases of the trajectory outlined, to arrive at a desirable 
altitude, velocity and attitude, is complicated because of the external forces 
Seine on the vehicle. The major forces affecting the vehicle during these 
phases are thrust, the earth's gravitational force, and the aerodynamic forces 
Of 1if% and drag, which act in a direction perpendicular and parallel to the 
instantaneous direction of flight, respectively. Gravitational and drag 
forces acting on the vehicle result in velocity losses during the flight and 
1965 detract from the efficiency of thé launch. The lift force may in certain 
cases be beneficial in that it may aid in turning the vehicle. 

For a specific vehicle, the important trajectory design parameters 
are velocity, altitude, and flight path angle at burnout. It is only possible 
however, to compute trajectory characteristics by numerical integration of the 


equations of motion from specified initial conditions. In this paper numerous 
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trajectories are developed, sine vehicle characteristics which are sappyox- 
imately representative of large chemical rockets of contemporary design, and 
varying the time of vertical flight and the maximum tilt angle. The effects 
of variations in two important vehicle design parameters are also included: 


namely, the initial thrust-to-weight ratio, n., and the mass ratio; the value 


јә 
of n; is introduced as an additional initial variable, while with the assump- 
tion of constant mass flow every point in each computed trajectory corresponds 
to burnout for some specific mass ratio. The burnout conditions are then 
examined as functions of the initial variables by using burnout angle as a 
governing parameter and cross-plotting. The nature of trajectory optimization 


to maximize burnout velocity or burnout energy is of particular concern to 


this study. 





CHAPTER 2 


VEHICLE DESCRIPTION AND AERODYNAMICS 


This study is intended to derive conclusions applicable to rocket 
vehicles similar to contemporary long range ballistic missiles, satellite 
launchers, and space vehicle boosters. Development of trajectory data re- 
quires the use of certain vehicle design parameters which identify aerody- 
namic and engine performance. To simplify preliminary work the "high-drag" 
configuration missile of Ref. 1 is selected as a model. It is believed 
that this design has aerodynamic characteristics representative of the class 
of vehicles described above. Rocket engine specific impulse is taken to be 
300 lb-sec/1b, and, again for simplification, this value is considered con- 
stant throughout the flight regime. A value of 10 is selected for the 
ratio of initial weight to burnout weight, defined as the mass ratio. This 
makes the propellant factor, the ratio of initial fuel weight to initial 
vehicle weight, equal to .9. Figure 1 shows the physical dimensions of the 
selected vehicle. 

Since the vehicle of this study is similar to the high drag missile 
of Ref. 1, the aerodynamic coefficients utilized are extracted from this 
source. A detailed explanation of the methods and procedures used to arrive 
at these values are set forth in Appendix G of that report. 

The missile, being axially-symmetric, has only a drag force imposed 


on it during the vertical flight phase and the gravity turn phase, since 





2 
W,/A = 3,010 1b/£t 


S/A = 11.0 


Fig. 1 Missile configuration. 


236,700 lbs 


300 sec 





(Reference 1) 
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during these phases the angle-of -acvcack 5 zero. her zero=) 1 i lvarap. [one 
15 made up Of three paris  baserdrap, skin friction: drag. and Tormrde>r: 
Mee ZCrO=115U drag coefficient; based upon both theoretical and enpiricar 
data, and representing the sum of these forces, is plotted versus Mach num- 
enfin Ејр, О. 

In order to simplify computer: programming, Һе curvevor D versus 
hasmenumber is divided into five segments. Each segment of the curve is 
tien represented by a straight line function. Тһе breakdown of the De 
curve and the approximating straight line functions are shown in Table 1. 
This straight line approximation of the curve representing E while being 
an approximation, is considered to be sufficiently accurate for the problem 
at hand. 

During the tilting phase of the trajectory, the missile is sub- 
5 Eko lift forces, as well as drag forces, since the missile has an angie 
of attack during transition from the vertical flight phase to the zgero-lift 
phase. Reference 1 outlines a cross-flow method of predicting lift and drag 
on bodies of revolution at an angle of attack. In this method the flow over 
the missile is separated into two components: one along the axial direction 
of the body, and one component normal to the axis. The axial flow exerts a 
zorceson the body in the axial direction while the cross flow exerts a force 
il anormal direction. Reference l derives equations for Cr, and Cp using 
this theory of cross and axial flow. 
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TABLA 1 


STRAIGHT LINE APPROXIMATIONS OF THE ZERO-LIFT DRAG COEFFICIENT CURVE 


FOR VARIOUS VALUES OF MACH NUMBE 





Do 
- len 
üzə 
0.72 
850M - .482 
1.0 
.983 - .323M 
‚522 - ,OĞOM 
‚328 - .023M 2 
5177 
TABLE 11 


STRAIGHT LINE APPROXIMATIONS OF THE CROSS-I'LOW DRAG CORFFICIENT CURVE 


FOR VARIOUS VALUES OF MACH NUMBER 


en wae —— 


M E 
Do 
О 
dy 
3.23M, ~ ‚49 
= 


ze an 


A arr 
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The term . is a drag coefficient due to the cross flow component. A plot 
of ... versus Mach number, assumed to appiy for this study, is shovn in 

Fig. 3. This plot is a series of straight line approximations of the cross- 
flow drag characteristics derived in Ref. 1 for the missile configuration of 
Fig. 1. These straight line approximations are described by functions as 
set forth in Table If. The straight line approximations are considered 
Subraciently accurate since, as it can be seen from the above equation, the 


erect of Cy is minor at small angles of attack. 
e 
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CHAPTER 3 


TRAJECTORY ANALYSIS 


The portion of the powered flight trajectory of interest in this 
study is considered in three phases, as discussed in Chapter 1. The first 
phase or vertical flight regime is followed by the tilt phase during which 
the vehicle is tilted from the vertical at a rate of two degrees per second. 

In Ref. 1 the tilt phase is approximated by impulsive tilting to 
5.5 degrees from the vertical during a one second time interval at the end 
of vertical flight time, followed by a gravity turn which continues until 
the desired conditions of attitude, altitude, and velocity are reached. The 
vehicle is assumed to be in the gravity turn as soon as the impulsive tilt- 
ing is accomplished. The impulsive tilting during a one second time interval 
Momjusvitied by determining that the required vehicle response time is less 
55 none second for the 5.5 degree tilt angle. This computation is made on 
the basis of the time required to tilt through the specified angle vith the 
maximum tilting moment available acting on the moment of inertia of the 
vehicle. For the present study, which is concerned with higher values of n; 
and consequent higher dynamic pressures during tilting, a tilt rate of 2 
degrees per second is selected as a reasonable maximum value. This is per- 
haps lower than necessary for tilting at sea level, but to have a basis for 
comparison, this rate is used for all trajectories computed. In the computer 


chis tilt rate is approximated by increasing the tilt angle, U, two degrees 





15 


per second until the tilt angle reaches the specified maximum programmed 


tilt angle, U,. This value of tilt angle is then held constant until the 
angle-of-attack of the missile becomes zero. At this time the tilt phase 
ends and the zero angle-of-attack or gravity turn phase begins. During 
the tilt phase thrust is considered to act parallel to the vehicle axis. 
The component of thrust required for tilt is considered a negligible loss 
compared to the total thrust vector. 

In the zero angle-of-attack phase thrust acts in the direction of 
£ inetantaneous velocity vector, vhich is also parallel to the missile 
axis. Turning is accomplished by the action of the earth's gravitational 
field. This part of the trajectory would logically be followed by a con- 
stant attitude or a "linear with time" thrust program, depending on the 
mission of the vehicle. In this study the gravity turn is continued until 
ninety percent of the missile mass is consumed. Since this paper deals only 


with single stage characteristics, staging is not considered and all results 


pertain to first-stage values. 
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CHAPTER İ 


EQUATIONS OF MOTION 


The equations of motion are developed using an inertial X, Y 
coordinate frame. This assumes a "flat", non-rotating earth, which is a 
good approximation during the early powered-flight phase of the type of 
7: .  vehicle considered. The gravitational acceleration due to the earth 
is assumed to be constant during the portion of the trajectory of interest 
in this paper. This also is a reasonable assumption when the altitude 
reached is small compared with the radius of the earth, as it is in this 
study. 

Rocket engine characteristics are simplified by assuming constant 
thrust and constant mass flow rate. Both of these quantities usually vary 
with atmospheric pressure, thrust increasing and mass flow rate decreasing 
as altitude is increased. This means that the specific impulse actually in- 
creases with altitude and that the initial thrust-to-weight ratio is based 
on the lower level of thrust found at sea level. The simplifications made 
in this study specify a constant specific impulse of 300 seconds, which may 
be thought of as representing an average value. The initial thrust-to-weight 
ratio in this study is therefore somewhat larger than it would be for an 


actual vehicle of comparable performance. 





wherein 


15 


Considering the vehicle as a point mass the equations of motion are 
D”Y/DT” » (F/m)cos U - gaye -(D/m)sind -(L/m)cos Y (1) 


D”X/DT” s (F/m)sin U -(D/m)cos $ 4- (L/m)sin)f (2) 


F/m = thrust per unit mass 


MA gravitational acceleration due to the earth 
D/m = drag per unit mass 
L/m = lift per unit mass 


The lift terms are considered positive in sign for the negative angle-of- 


attack condition which occurs during the tilt phase of the trajectory. Figure 


4 shows the vector relationships involved. 


Thrust, lift, and drag forces per umit mass are computed using the 


nomenclature of Ref. 2, in which T, is defined as a fictitious time when the 


total mass of the vehicle would be consumed. 


a w/m = 1./ni (3) 


From the conventional relationships between rocket parameters 


I, = F/w (1) 
m = W/g (1 - 1/7) (5) 
c z 1.80 (6) 


Applying (4), (5), and (6) to the various accelerations due to thrust, lift, 


and drag in (1) and (2) gives 


F/m = е/(Т, – Т) (7) 





16 





Fig. 4 Simplified diagram of vector quantities associated with 
the missile. 





D/m 


L/m 


= (SeVC Ag 7) /M, (Ty-T) 


= ($ CVC Ag 7, )/Wy (Dy-T) 
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(8) 


(9) 
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CHAPTER 5 


COMPUTER PROGRAM 


The computer used in this study is an IBM 650 digital computer 
located in the computation center of the Instrumentation Laboratory of the 
Massachusetts Institute of Technology. Although not comparable in speed to 
the larger digital computers such as the IPM 704, it is adequate for this 
study, computing an average trajectory in about ten minutes. The computer 
program is prepared using the MAC programming system developed by the 
Instrumentation Laboratory computation center. 

Time intervals for integration are varied according to the rhase 
255 7 creajectory. During the vertical flight phase the time interval is 
әсәб tour seconds for vertical flight times of four seconds and above, and 
emesoecond for vertical flight times less than four seconds. The time in- 
terval is reduced to one second during the tilt phase to maintain comparable 
accuracy in computing the rapidly changing trajectory quantities. At the 
completion of the tilt phase the time interval is increased again to four 
seconds and is held constant until burnout. 

The initial conditions for the equations of motion are set equal 
to zero for each run. Parameters held constent for all runs are: S/A, Ws /A, 
DU/DT, I,, and w. Variable parameters for each run are: T, T,, and Up. 


V 


IRo fietiticus burn-up time, T 


us equals I,/n;. Since I, is held constant, 


Ta is directly proportional to n Table III lists the numerical values of 


1 . 


the constants and parameters used. 
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TABLE III 


NUMERICAL VALUES OF CONSTANTS AND PARAMETERS 


Value Description 






Atmospheric density at sea 
level 


Gravitational conversion 
factor 


Gravitational acceleration 


acting on vehicle (assumed 
constant) 


1120 Ratio of planform area to cross 
section area 










Ratio of initial weight to 
cross section area 


2, 


300 sec 







Tilt rate 






Average specific impulse 
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Tilt rate is held constant at two degrees per second until Ue is 
reached. This is mechanized on the computer by increasing U instantaneously 
at the beginning of each one second time interval until U equals U At 
this point U, is held constant until the angle of attack becomes zero. Lift 
and drag are computed during the tilt phase in the manner shown in Chapters 
2 and 4. These calculations are made at the beginning of each time interval 
and are integrated as constants within the differential equation loop of the 
program. The error introduced by this approximation was small, as a result 
of the selection of integration intervals; in general the change in aerody- 
namic force from one interval to the next did not exceed three percent. 

When the angle-of-attack becomes zero, U is set equal to (90° - 7), and 
thereafter varies directly with 8 . This point marks the beginning of the 
gravity turn phase, 

In the gravity turn phase lift is set equal to zero and drag is 
calculated in the same manner as above using the zero angle-of-attack drag 
coefficient. The program is terminated when T equals .9 Ty which corres- 
ponds to the mass ratio of ten mentioned in Chapter 2. Values of velocity, 
altitude, range, flight path angle, and energy are punched for each computer 


time interval, and for burnout. 
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CHAPTER 6 


RESULTS 


For the study and results as presented here, all vehicle and 
trajectory parameters are held constant, except the initial thrust-to- 
werens ratio of the vehicle, the time of vertical flight, and the maxi- 
mum programmed tilt angle. The initial thrust-to-weight ratios used are 
15757 0, 2.5, and 3.0. The vertical flight times are varied from l 
second to 24 seconds, and the maximum programmed tilt angle is varied 
from -2 to 90 degrees. The trajectory calculations are continued in all 
cases for a total time, T = 0.9T,, i.e., a mass ratio of 10. 

Approximately one hundred trajectories were computed for this 
paper. Table IV lists values of mass ratio, velocity, altitude, flight 
ии ис, energy, scalar velocity loss due to drag, and scalar velocity 
ии Со gravity for some of the more useful trajectories. The effects 


of varying the parameters, n., T,, and U 


ol m, are shown in Figs. 5 through 18, 


which display burnout values of velocity, altitude, energy, and flight 
path angle plotted against maximum programmed tilt angle for a mass ratio 
of 10. Separate plots are shown for each vertical flight time used. 

Scalar velocity loss due to drag is shown in Figs.19 through 22. 
Vp, divided by the ballistic coefficient, ne is plotted versus burnout 
flight path angle for each n; and each Ty. A more general presentation is 
Gobainea with the ballistic coefficient, whieh uses the vaiue of “De at Mach 


240: 





— 


Since side loading is an important consideration in large rocket 
design, plots of the maximum lift load factor versus Ди for cach H and each 
T, are included as Figs. 23 through 26. 

The length of time that the rocket is subjected to lift loads is 
also of interest. Figure 27 shows the time required to tilt the missile so 
that it will attain a burnout flight path angle of thirty degrees. This is 
plotted against nį to show the large increase in time required for tilting 
аз п, is increased. 

An optimization study is made, based on finding the combination of 
parameters which would give the highest specific energy at burnout for each 
n, investigated. Energy is first maximized for fixed values of ð, by plot- 
ting energy versus the value of U. corresponding to particular vertical 
flight times. The maximum energy points are then cross-plotted against ð y, 
and the value of Un corresponding to the maximum energy point for the fixed 
values of S p- As a cross-check, this procedure is reversed so that energy 
is maximized for fixed values of Un and cross-plotted in the same manner. 
Also included in these figures are the values of burnout velocity and burn- 
out altitude which occur at the maximum energy points. The optimization 
results are shown in Figs. 29 through 31. 

Representative trajectories for various values of n, are identified 
in detail in Table V. These particular trajectories were chosen because they 


were closest to the maximum energy trajectories for each case. 
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16 
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18° 
L41 
IGL 
180 
147 
142 
162 
180 


17” 
L41 
161 


26" 
142 
162 
180 


a 
139 


139 
180 


29” 
Tı 
161 
180 


dı 
139 


152 
180 


TABLE IV 


TRAJECTORIES INVESTIGATED SHOWING BURNOUT 


VALUES FOR VARIOUS VALUES OF MASS RATIO 


1.09 
3.39 
2.13 


1.15 
3.45 
5.26 
LOMO 


12.16 
5 
4.88 
10.0 


ar 
3.39 
2.13 
10.0 


Tor 
3.28 
4.88 
TORO 


5 
» 232 
7,190 
10,684 
16,700 


252 

8,169 
12,086 
18,130 


314 
1,843 
2,970 


506 
7,20% 
10,731 
16, 400 


25 
275 
10,565 
17,100 
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Fig. 5 Variation of Vy with 7 for various values of ni at a 
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Fig. 14 Variation of E, with U, for various values of n, at a T, 
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Fig. 16 Variation of By, with U, for various values of n. at a 
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Fig. 17 Variation of 3, with ue for various values of ni at a 
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Fig. 18  Variation of L. with U„ for various values of n, at a T 
of 24 seconds. 
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Fig. 25 Variation of maximum lift load factor with E for various 
values of n, at a T, of 16 seconds. 
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under maximum energy conditions for n, of 2.0. 
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Fig. 29 Values of Un and Tr required to obtain a specified R. 


under maximum energy conditions for m Of “2.0. 





56 





io 
Ex 10. 
(ft-10/slug) 
1 
ik 
80 
60 
0 ДО 
m 
(deg) 
20 
0 
2} 
16 
T 
v 
(sec) 8 
O 
90 60 30 O 
ду (deg) 


Fig. 30 Values of U, and T, required to obtain a specified X 


under maximum energy conditions Tor H. 01557 is 





700 
600 
500 

BG 1073 

(ft) 100 


300 


200 


19 


18 


Vp X 10, 
(fps) 


11 


16 


90 60 30 
Y, (aeg) 


Tig? 01 Verves or Ya, and Vy at maximum energy. 





21 








* End of tilt phase. 





108 


* End of tilt phase. 


COMPUTER RESULTS FOR REPRESENTATIVE TRAJECTORIES 


265 . 78352 
546.99767 
1152.0408 
1227.7526 
1933.2139 
2087.3356 
2731.9199 
3626.1614 
4799.5585 
6325.6918 
8308.6788 


10972. 134- 


14886. 326 
18668 . 867 


199.15440 
409.79788 
632.93575 
870 . 24830 
992.92178 
1244. 3838 
1705 . 8846 
2151.1125 
2721.8130 
3476.8072 
4427. 3496 
5617.1983 
11.00.7933 
8964.8019 
11391.601 
14785 .568 
18585 . 407 


525.70761 
2120 7717 
8609.4335 
9726. 6649 
14867 .757 
28129.005 
45316.107 
67522.151 
96429.217 
134027.09 
182846. 45 
246406. 84 
330622. 30 
398806. 33 


ә“ 


393.99420 
1593.1159 
3600.1614 
6387.4618 
8071.9079 
12031.649 
2213.8303 
3464. 219 
49616.932 
0712 121 
89914.139 
116914.58 
14.9892 .68 
190254. 90 
240036.72 
302500. 64 
361322.79 


TABLE V 


T = 1 sec 


У 


deg 


87.002181 
80. 949374 
70.776125 
70.160493 
69.635561 
66.455353 
63.757048 
61.496942 
59.667247 
58.208353 
512055556 
56.151720 
55.454918 
54.922312 


T, = 1 sec 


y 


deg 


86.681749 
80.054237 
73.392714 
66.923305 
64.713429 
63.844218 
58.869145 
54.543873 
Do r 
47.445993 
4l.. 721666 
42.479115 
40 . 646909 
S RTE 
37.954541 
36.996412 
36.152776 


- O 
Ün - 20 


E 


ft-1b/slug 


52234.583 
218126.61 
940599. 34 
1066634. 4 
1656889.0 
3083509.0 
5189695 .9 
8746984.1 
14620399.0 
24319383.0 
HO 399983.0 
68121766.0 
121438810 
187094520 


Ј = 269 


m 


E 


ft-1b/slug 


32507.627 
135225.11 
316135 .60 
584176.57 
752652.80 
1161352. 4 
2167300.1 
3428287. 3 
5300510.7 
8224919. 4 
12693614.0 
19538075 .0 
30033288.0 
46305429.0 
72607238.0 
119039180 
134333900 
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Ур 
fps 


‚20644683 
IE 
27.850098 
37. 685454 
87.560342 
331. 30749 
610.93239 
802. 44887 
923 2300 
982.43218 
1000.4780 
1005 .5098 
1006. 2890 
1006. 2980 


Vp 


fps 


.11541601 
1.227409 
4 .6832111 
11.738844 
18.598280 
39.663503 
180.90261 
425 .59583 
647. 48641 
808.12201 
924.. 52635 
994. 29298 
1026.7656 
1040. 6460 
1045. 6924 
1047 .0206 
1047 .0806 





TABLE V (Continued) 


= = o 
© mo T = 16 вес Un 18 

T V Y y E Vp 
sec fps et deg ft-1b/slug fps 

Y 67.000877 132.68853 89.999999 6513.6862 

8 137.96242 541.24719 89.999999 26930 .929 .060251961 
12 212.91698 1241.5799 89.999999 62613.474 . 31773094 
16 291.88676 2249.6993 89.999999 114980.88 . 93047546 
20 374.09575 3580.1953 88.110768 185163. 20 2. 3827188 
24 459.54522 5236.3599 82.724083 274065 .82 4 .8911670 
28 551.81331 7208.5201 76.258076 384176.26 8.292982 
” 626.59842 8900.2976 72.774338 482671.41 11:20 07 
35 732.62606 111454.434 71.827600 636906 . 00 15.982745 
43 964. 12244 17603. 648 65.552245 1031146.7 31.150525 
51 1207. 2451 25045. 312 59.270032 1534529.5 72. 764797 
59 1450 ,8226 33523.212 5 3.090988 2131020 .6 157.13915 
67 1706.2258 42720.706 17.055876 2830101.4 278.08037 
75 2027.7450 52525.718 41.291939 3745840.0 386.30959 
83 2423.2376 62935 . 490 35.985004 4960929.9 478.66841 
91 2898.h151 73921.632 31.228459 6578763. 4 554.25915 
99 3452.0813 85422.589 27 .038602 8706823.1 619.84834 
107 4089.9748 97389.100 23. 380049 11496062.0 676.75458 
115 1821.4522 109613.97 20.200540 15149926.0 724.96740 
123 5658.5944 122137.53 17.14 44041 19939505 .0 765.06602 
131 6621.16994 134856.15 15.057469 26262321.0 794. 20610 
135 7154.6835 1171776 13.987885 30140033.0 807 . 78328 
147 9009. 8838 160685. 28 11.208726 45758899.0 844. 95417 
15 10523.848 173759.37 9.6691722 60066233.0 867.32021 
163 12354.673 186986. 81 8. 3463049 82335093.0 888.65247 
171 14652.533 2004.90.62 72179499 11379896 .0 910.07619 
180 15181.017 216305 .47 5.8239900 17223411.0 936.87476 


* End of tilt phase. 
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* End of tilt phase. 


fps 


132. 98571 
274 .54739 
125 .70506 
755.87968 
1118.9996 
1473.1912 
17.97.0701 
2184.0412 
2687 . 623} 
0 12.120L 
4064 .6675 
4967.0819 
6045.2513 
7236.8161 
8894.5880 
10821. 625 
13307. 439 
16759.873 
18184 .668 
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TABLE V (Continued) 


2 | 12 

m, s ar Tara 16 sec u: 
Y у E Vp 
It deg ft-1b/slug fps 

262.85677 86.04007 17299.766 .051083963 
1061.2315 18.231773 71832.250 .53748826 
2406.9641 76.811880 168054.18 1.5548h19 
6784. 949 69.191874 503976.21 8.4062369 
13404.599 62.519935 1057360. 3 31.113584 
22048.882 56.628937 1794548 .0 117.69074 
32228 . 993 51.247801 2653105 .7 295.59529 
43659. 480 46.242080 3789720.3 481.84897 
56560.978 41.736274 5431455 .5 629. 68643 
“20 9,7 37.809023 TT16023.5 746.65796 
87822.151 34.441157 11086355 .0 840.03946 
106668.91 31.587084 15767923.0 907. 74155 
128059. 42 29.174504 22392722.0 951.86607 
152389. 86 27.1423 31817434.0 976.00713 
180171.00 25.440375 15353679.0 991.09169 
212127.15 24.016065 65378775.0 999. 40505 
249378. 99 22.835027 96567497.0 1003.5360 
293825.31 21.546874 149900220. 1005 .0832 
312997. 36 21.250934 180911870. 1005 .2250 





CHAPTER T 


DISCUSSION OF RESULTS 


The status of the vehicle at burnout is of prime importance. 

This vehicle status is best described by the burnout quantities, Vp? Ip Ep» 
and § 1, as displayed in Figs. 5 through 18. 

These burnout values are not shown for the vehicle where the value 
of H is 1.5 and the T,, is 1 second, due to the fact that the vehicle passes 
through the horizontal and heads back towards the earth before reaching 
burnout. This result is readily explained since the vehicle is relatively 
slow and turns at a low altitude where drag losses are very large. 

The vehicle does reach burnout conditions for the other values of 


2777 7 er, For the case vhere n, is 1.5, the value of Va, peaks at a very 


V? al 
low value of U The larger values of um cause the relatively slow vehicle 
to turn more at a low altitude. This reduces Ҹу due to the fact that the 
vehicle operates for a longer period in dense atmosphere where the drag 
velocity loss is very large. 

As In is increased, the maximum value of Vy, obtained for an n4 of 
1.5 Occurs at the higher values of Um: This happens as a result of the in- 
ereasea velocity and the higher altitude reached before the tilting is 
commenced. Drag velocity loss is much less under these circumstances. It 


is true that gravity velocity loss increases as T; is increased but it does 


not offset the reduction of the velocity loss due to drag. This fact is 





further borne out by observing that the curves indicate higher burnout alti- 
cudes are reached as the value of un is reduced. 

The burnout altitude reached varies inversely with the value of U, 
and directly with the value of r for the vehicle with an ni OBTL: BODA 
of these phenomena are readily explained since the altitude attained is a 
direct function of the vertical component of the velocity vector which is 
directly affected by these two factors. 

The value of O b of the vehicle for the values of n; ~ 1.5 generally 
decreases with increased values of U, and with an increase in Т,. A large 
U,, allows a greater amount of turn of the vehicle, hence a smaller ¥ 4: At 
the higher values of T, the vehicle has a greater velocity before commencing 
the turn, and hence does not get turned as much since the gravity vector 
causing the turn after um is reached is small in comparison to the vertical 
component of the velocity vector of the vehicle. 

Generally, the burnout energy of the vehicle follows the trend of 
the velocity at burnout. This is explained by the fact that energy is directly 
proportional to the square of the velocity. 

As the value of n, is increased the burnout velocity attained be- 
comes less and less dependent upon the value of Uns which is indicated by the 
fact that the curves of Vy versus Us tend to flatten out as the value of ni 
increases. Vertical flight time does not affect the burnout velocity to any 
great degree and at higher vertical flight times the value of U. appears to 
have less affect on the L reached. The basic reason for these effects is 
that with a high ni and T,, the vehicle is out of the very dense atmosphere 
before turning, and hence the drag velocity loss is relatively small. 

As H, is increased the burnout altitudes become less and less at 


the low values of U, and just slightly more at the high values of Un At 
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the higher vaiues of О the high H, vehicle does not get turned as much as 
the lower n, vehicle and hence attains a Slightly higher burnout altitude. 

Vertical flight time seems to have little effect on the burnout 
velocity attained. It appears to be a slight factor at low values of Т,, 
but for values of T, greater than 8 seconds, vertical flight time has no 
appreciable effect upon the burnout velocity reached. The value of бу де- 
creases with an increase in ni for the reasons previously explained. 

The drag velocity loss of the missile for the various trajectories 
stuđied is shown in Figs. 24 through 27. The variation of this loss with 6 T 
and n; for four values of 1. is shown. The more prominent indications of 
these results are: 
(1) High n; missiles have the highest drag velocity loss. 
(2) The drag velocity losses decrease slightly with an increase 

in vervical flight time. 


(3) Drag velocity losses greatly increase as the missile attitude 


approaches the horizontal at burnout. 


The higher drag velocity losses accompanying an increase in the 
value of n, 1S a result of the hiecher7velscitiessartained by Lhe vehicle aL 
Over altitudes. Since drag force is directly proportional to the square of 
the velocity and to the density of the atmosphere, this force, and hence the 
resulting velocity loss, is large for a high n; missile. Since the drag 


velocity loss may be expressed as: 


= D 
Vp = g, f? at 


it must be realized that a slow, large vehicle would experience a lower drag 


velocity loss than would a smaller, faster missile. 


m CY 0 





h 


The drag velocity loss decreases slightly with an increase in ver- 
tical flight time because the missile is at a higher altitude and hence in 
a Less dense atmosphere before it commences to turn. It therefore spends 
less time in the denser atmosphere of low altitudes. This decrease in drag 
velocity loss as indicated in the aforementioned plots, is not as great as 
would be expected. The drag loss is plotted versus Y b), and to get to the 
same value of gy for a high vertical flight time as for a low vertical 
flight time, the vehicle must be turning at an angle of attack for a longer 
period of time since the missile has a greater velocity at the start of the 
tilting phase. This increased time of tilt with an angle of attack increases 
the drag coefficient due to the induced drag present while this condition 
exists. 

Nearly the same reasoning applies to the condition of increased 
drag velocity loss for a smaller value of Bu attained for a particular 
vehicle with a given n, at a certain Т,. The smaller value of LA requires 
that the missile be tilted at an angle of attack for a longer period and 
hence the drag coefficient is again increased. 

The values of ny, which are plotted against U, in Figs. 28 through 
31 show the maximum negative lift load factor encountered during a trajectory 
of specified T, and Un: Lift is always negative in the tilt phase, that is, 
in the direction of rotation of the vehicle. Lift load factor increases 
rapidly with increased Ye for Lie snarh ni trajectories. The slope of the 
lift load factor curve increases as Т, 159 increased, as could be expected 
from the higher dynamic pressures caused by the higher velocities associated 
with long vertical flight times. An interesting aspect of the higher İV 
plovs is that the maximum Hr. seems to level off, reaching a maximum of about 


1.2 g's for an ni of 3.0. In these cases the missile has reached the less 





ни 
e 
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dense portions of the atmosphere, where the extremely low density has offset 
the increased velocity, and the maximum ny encountered has become essentially 
constant. The large values of ny encountered for high nį trajectories and 
the accompanying bending moments are too great for most contemporary liquid 
fueled vehicles. Use of trajectories of this nature require the heavier 
structural design associated with solid fueled vehicles. 

It might be mentioned at this point that the rate of tilt should 
have a considerable affect on lift. This program uses a tilting rate of two 
degrees per second, which is considered a nominal rate for a large rocket- 
powered vehicle control system to achieve, but a minimum rate to reach the 
maximum programmed tilt angle within a reasonable time. The two degree per 
second tilt rate was selected to keep the lift load factor within practical 
bounds. A study of methods for obtaining minimum lift loads through differ- 
ene types of tilting programs is of interest, but beyond the purposes of 
this paper. 

The variation of tilt time necessary to reach a 8, of 309° for 
four different vertical flight times is plotted against n,; in Fig. 27. The 
Curves indicate that the tilt time increases as both Т, and n, increases. 

The higher the value of T, used, the greater the velocity of the 
vehicle becomes before the vehicle starts to turn. Since this is the case, 
it takes much longer to turn the vehicle from the vertical position to a Ús 
of 30° at the higher value of T; than at the lower values. This same type 
of reasoning may also be applied to the second observation, in that at high 
values of ni the vehicle gains greater velocities sooner than at low values 
of H It therefore requires a greater tilt time to reach a 8, of 30° at 


the higher values of e 





66 


The optimization program undertaken in this paper is based on max- 
imum specific energy. The combination of Т, әла Un which would give maximum 
burnout energy for a specified burnout angle are determined. Energy is max- 
imized in this manner for the three higher values of N; The resulting values 
of Ey, Un and Ty are plotted versus burnout angle in Figs. 19 through 2l; 

t is not possible to obtain a good optimization for an H, of 1.5 witha the 
data available. At this value of n;, burnout energy is extremely sensitive 
to Т, ала У,, and any attempt at optimization requires a large number of tra- 
jectories. Accordingly, no optimization is included for an ni OL Lis 

Good maximum energy points are found for burnout angles below about 
thirty degrees. In the range of thirty to fifty degrees, two approximately 
equal maximum energy points appear. One of these points occurs at a value 
of A, consistent with the maximum points found in the thirty degree and below 
range, while the other occurs at the minimum Т,, which is one second. Above 
the thirty to fifty degree range the maximum energy points are found at the 
minimum values of T and U, for the particular burnout angle. These same 
characteristics, in varying degree, are found for each value of Nj- The 
overlapping in the plots of T, and Un versus burnout angle in the thirty to 
fifty degree range show that in this area maximum energy can be obtained by 
using either of two combinations of T, and Um: This effect is undoubtedly 
caused by the non-linear action of drag on the trajectories. A point is reached 
for each ni where the beneficial effects of early tilting, and consequent earl- 
ier alignment of thrust and velocity vectors in the desired direction, are 
overcome by the higher drag losses associated with large programmed tilt angles 
at low altitudes. At this point it is necessary to use a period of vertical 
flight time to get the vehicle out of the denser portions of the atmosphere 


during the tilting maneuver. 
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The amount of vertical flight time required to maximize burnout 
energy for low burnout angles does not increase indefinitely as burnout angle 
approaches zero, but tends to level off in the neighborhood of 16 to 20 seconds. 
This is especially evident in Fig. 21 where nj 1s 3.0.9) the High velocity 
reached at the end of a long vertical flight time causes an increased negative 
lift load during the tilting phase, which increases the drag coefficient, 
“hereby reducing burnout velocity and burnout energy. 

It is interesting to note that at burnout angles of about ten degrees 
and beiow, the maximum burnout energies for all three values of m fall very 
close together. This indicates that for very low burnout angles the advantage 
of a high initial acceleration rocket is questionable, since the same burnout 
energy can be obtained using a lower initial acceleration with lower aerodynamic 
loads. At other values of burnout angle the higher burnout energies obtained 
from high m rockets are apparent. The point at which maximum overall burnout 
energy is reached starts at a burnout angle in the vicinity of twenty-five 
dor eg for an n, of 2.0, and moves in the direction of increasing burnout 
angle as ni increases. In this case, the increased drag associated with high 
initial acceleration and low burnout angle causes the maximum overall energy 
points to fall at higher burnout angles for the higher values of n,. 

Values of burnout velocity and altitude for the maximum energy condi- 
items Giscussed before are shown in Fig. 22. It logically follows that, since 
energy is a function of velocity and altitude, these curves are coincident in 
the same manner as the maximum energy curves at low burncut angles. At higher 
burnout angles, burnout velocity increases and burnout altitude decreases as 
ns 1s increased., 

Each representative trajectory tabulated in Table V is the nearest 


trajectory available to the overall maximum energy case for that particular nj: 
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This is the optimum trajectory for the burnout angle listed only, and it can- 
not be said that it is the optimum trajectory for any other attitude reached 
before burnout. Naturally it is not applicable to a burnout angle lower than 
that listed. The question arises whether or not more energy is obtained by 
using the tabulated trajectory for the overall maximum energy case until the 
desired attitude is obtained, followed by constant attitude thrust until 
burnout, than by using the maximum energy trajectory for the burnout angle 
corresponding to the attitude desired. It seems that if the vehicle is above 
“һе denser atmosphere, the energy generated after constant attitude thrust 

is started would be about the same as that for the gravity turn. A constant 
attitude thrust program would give higher altitude with lower velocity than 
the gravity turn, although the actual difference between the two programs 
would depend upon the time of application of the constant attitude thrust 
program. If the vehicle reaches the desired attitude in the early tilt phase, 
before leaving the denser atmosphere, a constant attitude thrust program 
would involve lift loads and additional drag, but it would get the vehicle 
out of the sensible atmosphere sooner. In either case the non-linearity of 


the problem would require a separate computation for each situation. 
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CHAPTER 9 


CONCLUSIONS 


Determining the optimum powered flight trajectory for a large 
rocket powered vehicle is a complex problem, which is strongly influenced 
by the desired trajectory burnout angle and the rocket initial thrust-to- 
weight ratio. The burnout angle may be considered a design parameter for 
the booster trajectory, since different burnout angles are required for dif- 
ferent missions. This paper shows that a combination of vertical flight 
time and initial tilt angle to give maximum energy at burnout can be deter- 
mad forany desired burnout angle. There is, in addition, one value of 
burnout angle, which gives maximum burnout energy, for each value of initial 
acceleration. In this manner an optimum booster flight trajectory is avail- 
able for any'desired burnout angle. 

Usually the vehicle with the higher initial acceleration will have 
the greater burnout energy. At low burnout angles, however, in the zero to 
ten degree range, values of burnout energy for a wide range of initial accel- 
erasions closely coincide. For low burnout angles, therefore, the initial 
acceleration of the vehicle is of little consequence with respect to maximum 
energy cptimization. In fact, it can be said that a lower initial accelera- 
tion is preferable for low burnout angles, since the high lift load factors 


associaved with high initial accelerations are avoided. 
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The time required to tilt the vehicle from the vertical to a point 
where the angle-of-attack is zero is also quite high for high acceleration 
vehicles. The procedure involving a relatively small initial tilt angle 
followed by a gravity turn to the desired burnout angle is no longer feasible 
with high initial accelerations, which leads to the conclusion that the tilt 
phase for high initial acceleration vehicles must be programmed in its 


entirety. 
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APPENDIX A 


ATMOSPHERIC DATA 


The atmospheric data used is based on the ARDC model atmosphere 
of 1959 as described in Ref. 3. In order to facilitate computer procedures 
the density ratio (€/@,) and sonic speed data are treated in a simplified 
manner. 

Utilizing the atmospheric density data of Ref. 3 a plot of the 
ratio of e@ /e, is made extending from sea level to an altitude of 400,000 
feet, as shown in Fig. 32. The resulting curve is divided into four seg- 
ments which are accurately approximated by appropriate exponential functions. 
The resulting segments of the curve and respective describing exponential 
functions representing the density ratios are shown in Table VI. 

sonic speed is plotted versus altitude from sea level to an alti- 
tude of 400,000 feet as shown in Fig. 33. The curve results in 8 series of 
five straight line segments. Straight line functions are used to describe 
De segments of the curve. Table VI displays these functions and their 


respective areas of applicability. 
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TABLE VI 


ATMOSPHERIC DATA APPROXIMATION FORMULAS 












Altitude Density Ratio Speed of Sound 
(ft) ES (fps) 
S Y/32,000 
e-Y/32, 1120-.00417 Y 
36,800 





1.65 e”Y/ 20,800 
82,500 
ҝӧ | 
168,000 
263,000 






813.78 1001077 








ти 846.5 
450,000 





300 


250 


200 


“x 108 
(ft) 


150 


100 


2. 
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Variation of speed of sound with altitude. 
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